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Ctase propose (M£^cL/td/) 
LMASSUinO 

Un metodo automatico di localizzazione ed estrazione del contorno di 
strutture present! in una immagine che usa una nuova classe di filtri, 
tiLerivati dal inomento centrale assoluto del primo ordine. II metodo prevede 
il filtraggio di una immagine di partenza tramite il momento centrale 
assoluto e(n,m) dell' intensita dei pixel di detta immagine, essendo n e in 
le coordinate di ciascun pixel, in cui il momento centrale assoluto e 
ottenuto con le seguenti fasi: determinazione per ogni n,m della media 
TLocale calcolata su un intorno del pixel di coordinate n,m dell' immagine 
d^ partenza, ottenendo una prima immagine filtrata; determinazione per 
-dl^x n,m della sommatoria delle differenze in valore assoluto era 
J^intensita del pixel di coordinate n,m della prima immagine filtrata e- 
JL'intensita di tutti i pixel contenuti in un intorno del pixel di 
tioordinate n,m di detta immagine di partenza, o di una seconda immagine 
jfiltrata da essa derivata. Le immagini da trattare possono essere di 
Trario tipo e provenienti da vari settori, tra cui robotica, concrollo di 
processi industrial!, medicina, applicazioni multimediali , sis^emi di 
3icurezza, e possono essere immagini a colori o a livelli di grigio. Le 
immagini bidimensionali possono dare luogo a immagini volumetriche o 
immagini panoramiche se acquisite come sequenze spaziali di sezioni 
bidimensionali. 
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Descrizione dell' invenzione industriale dal titolo: 
''METODO E APPARATO PER L'ESTRAZIONE DEI CONTORNI DELLE 
STRUTTURE PRESENTI IN UNA IMMAGINE MEDIANTE L'IMPIEGO DI 
UNA CLASSE DI FILTRI NON LINEARI" a noma di CNR CONSIGLIO 
NAZIONALE DELLE RICERCHE, di nazionalita italiana e con 
sede a Roma, e a . nome di ESAOTE S.p.A., di nazionalita 
italiana e con sede legale a Casale Monferrato (AL) . 

Invenrori designati: DEMI Marcello; GEMIGNANI 
Vincenzo; PATERNI Marco; BENASSI Antonio. 

DESCRIZIONE 
Ambiro dell^ invenzione 

La presente invenzione riguarda un metodo di 
localizzazicne ed estrazione dei contorni delle strutture 
presenti in una immagine. 

Inoi-re, 1' invenzione si riferisce ad un apparato in 
cui risiedono filtri non lineari applicabili in procedure 
di localizzazione ed estrazione di detti contorni. 

Le- irrmagini da trattare possono essere di vario tipo e 
provenienti da vari settori, tra cui robotica, controllo di 
processi industriali, medicina, applicazioni multimediali, 
sistemi di sicurezza, e possono essere immagini a colori o a 
livelli di grigio. 

Nel caso di applicazioni in medicina, in tal caso 

riguardanti 1' imaging biomedico, esse si riferisco 

!ng, Marco Celestino 
ABM Ag^nzia BreylieUi & Marchi 
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normalmente ad un organo e sono ottenute con varie tecniche 
come ad esempio impulsi ultrasonici, PET, SPECT, CAT, MR, 
etc. In particolare, possono essere immagini anatomiche, o 
immagini di funzione, ottenute per mezzo di sequenze 

5 temporal! di viste anatomiche di una zona particolare di un 
organo, o immagini di perfusione, ottenute sullo stesso 
organo dopo trattamento del paziente con sostanze che 
mettano in risalto la perfusione nell' organo. Oppure, 
possono essere immagini di tracciati (es. ECG) acquisite con 

10 uno scanner con lo scopo di convertire il tracciato su carta 
in segnale digitale. 

Le immagini bidimensionali possono dare luogo a 
immaqini volumetriche o immagini panoramiche se acquisite 
come sequenze spaziali di sezioni bidimensionali. 

15 Descrizione della tecnica nota 

Le procedure di estrazione dei contorni 
costituiscono un importante aspetto dell' elaborazione 
delle immagini, poiche sono correntemente utilizzate nei 
settori suddetti. II contorno, definite come una rapida 

20 variazione di intensita e/o colore tra due regioni 
adiacenti, pud essere estratto attraverso numerosi 
approcci tra cui quelli piu usati si basano su filtri 
derivativi di primo e secondo ordine. 

Filtri derivativi di primo ordine 
25 I filtri derivativi di primo ordine misurano la 

Ing, Marco Celesiino 
ABM Agenzia Breveiti <Sc Marchi 
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variazione dell' intensita luitiinosa in una predefinita area 
di interesse, Le discontinuita in questo caso sono 



filtro derivative. 

Uno dei piu usati e il Gradiente di Gaussiana (GoG) , 
la cui risposta e soddisf acente su contorni rettilinei privi 



in presenza di contorni con valori elevati di curvatura (per 
es. spigoli) , alle estremita di linee o sulle barre, alle 
intersezioni tra differenti discontinuita ed in punti 
isolati (spot). In presenza di rumore, I'uso di una soglia 
predefinita permette di discriminare i massimi locali spuri 
da quelli generati dalle discontinuita. Molte procedure 
utilizzano un unico e predefinito valore di soglia su tutta 
I'immagine. Questo tipo di approccio risulta tuttavia poco 
efficace poiche le caratteristiche dell' immagine 
(luminosita, contrasto, rumore, ecc.) non sono, in generale, 
costanti sull'intera area. Per esempio, le parti piu scure 
di una immagine sono spesso caratterizzate da variazioni di 
intensita minori rispetto alle aree piu luminose. Quindi, un 
valore di soglia adatto per le zone luminose pud non essere 
adatto per le zone piu scure e viceversa. Questi problemi 
hanno portato alio sviluppo di procedure di tecnica nota 
capaci di estrarre localmente i valori di soglia che meglio 
si adattano ad una specif ica area dell' immagine, associanc^d 



localizzate cercando i massimi locali della risposta del 



di incroci con altri contorni, ma tende a perdere efficacia 
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ad ogni pixel dell' iinmagine una soglia legata alle 
variazioni di luxninosita. In particolare: 

- US6094508 utilizza una soglia ricavata attraverso un 
modello che approssima e misura la risposta della visione 

5 umana alle variazioni di intensita locali; 

- EP0414505 utilizza una. soglia dinamica ricavata dalla 
stima dell'errore di diffusione nell' intorno di ciascun 
punt o dell ' immagine . 

- US6061476 utilizza un approccio nel quale, attraverso una 
10 tecnica di sottrazione, si cerca di ottenere una soglia che 

e proporzionale al contrasto tra il fondo e la 
caratteristica di interesse delle quali si desidera estrarre 
il contorno . 

- US4220972 risolve il problema di estrarre contorni d^ 
15 immagini con basso contrasto, utilizzando una sogltfe"^ 

proporzionale alia misura della luminosita locale^ 
dell' immagine. 

Filtri derivativi del secondo ordine. 

I filtri derivativi del secondo ordine producono uno 
20 zero-crossing in corrispondenza dei contorni delle strutture 
presenti nelle immagini. Questo zero-crossing e dato dal 
fatto che i contorni sono individuati come delle 
discontinuita, in corrispondenza delle quali 1' intensita 
della risposta attra versa il valore zero. 

25 11 filtro derivative del secondo ordine piu noto i 

fng. Marco Celestino 
ABM Agenzia Brevet ti & Marchi 
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letteratura e il Laplaciano di Gaussiana (LoG) , ampiamente 
utilizzato in applicazioni di computer vision. Esso e spesso 
approssimato con una Differenza di due Gaussiane (DoG) aventi 
differenti aperture. In questo caso la localizzazione dei 

5 contorni si basa sulla ricerca dei punti di zero-crossing 
della risposta del filtro derivative. Tale operazione non e 
facile, soprattutto in presenza di rumore, perche 
1' operazione derivata di secondo ordine esalta il rumore. 

Anche in tal caso, il rumore viene affrontato 

10 utilizzando un' operazione di sogliatura, come ad esempio in 
US5086219, adottando una soglia legata alio strenght dello 
zero crossing, ossia alia pendenza della curva nel punto di 
attraversamento dello zero; il presupposto e che, in 
generale, alle discontinuita corrispondono zero-crossing piu 

15 intensi (cioe attraversamenti dello zero con pendenze piu 
elevate) di quelli normalmente generati dal rumore, che in 
questo modo pud essere filtrate. 

In alternativa, US5142592 propone una procedura basata 
sul filtro LoG per la localizzazione di contorni vicini e 

20 paralleli; in questo caso, il vincolo stesso di parallelismo 
viene sfruttato per contrastare i problemi derivanti dalla 
presenza di rumore. 

In generale, i filtri derivativi di I e II ordine 
danno risposte piu scadenti in corrispondenza di angoli, 

25 spigoli, vertici, intersezioni, variazioni di curvatura^i^?^ 

Ing. Marco Celestino ' ^ 
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ecc. Cio comporta, oltre ad un risultato non del tutto 
soddisf acente, un aggravio computazionale . 

Sintesi dell ' invenzione 

6 scopo della presente invenzione fornire un metodo 
per la localizzazione e I'estrazione di contorni di 
strutture presenti in una immagine che consenta di ottenere 
una immagine con caratteristiche in generale superiori 
rispetto a quelle ottenute con i metodi di tecnica nota. 

un altro scopo della presente invenzione fornire 
un metodo per la localizzazione e I'estrazione di contorni 
di strutture presenti in una immagine che consenta di 
mettere in evidenza con elevata risoluzione punti chiave 
di un' immagine, quali angoli, spigoli, vertici, 
intersezioni, variazioni di curvatura, linee, punti 
isolati, ecc. 

E un ulteriore scopo della presente invenzione 
fornire un metodo per la localizzazione e I'estrazione di 
contorni di strutture presenti in una immagine che 
permetta di realizzare le operazioni di calcolo piu 
velocemente che con i sistemi di tecnica nota. 

E uno scopo particolare della presente invenzione 
fornire un metodo per la localizzazione e I'estrazione di 
contorni che permetta di distinguere il bordo interne dal 
bordo .esterno di strutture presenti in una immagine 
corrispdndenti a detti contorni. 



Ing. Marco Ceiestino 
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t anche scopo della presents invenzione introdurre una 
nuova classe di filtri, applicabili in tale metodo di 
localizzazione ed estrazione di contorni. 

L' invenzione verra ora sintetizzata facendo 
riferimento ad immagini mappate a livelli di grigio, 
essendo chiaro che essa puo essere estesa anche a immagini 
a colori. In particolare, si parlera di intensita sia per 
definire 1' intensita dei livelli di grigio, che dei colori. 

Secondo 1' invenzione, i suddetti scopi sono 
raggiunti da un metodo automatico di localizzazione ed 
estrazione del contorno di strutture presenti in una 
immagine che usa una nuova classe di filtri, derivati dal 
momento centrale assoluto del primo ordine. 

II metodo prevede il filtraggio di una immagine di 
partenza tramite il momento centrale assoluto e(n,m) 
dell' intensita dei pixel di detta immagine, essendo n e m 
le coordinate di ciascun pixel, in cui il momento centrale 
assoluto e ottenuto con le seguenti fasi: 

- determinazione per ogni n,m della media locale 
calcolata su un intorno del pixel di coordinate n,m 
dell' immagine di partenza, ottenendo una prima immagine 
f iltrata; 

- determinazione per ogni n,m della sommatoria delle 
differenze in valore assoluto tra 1' intensita del pixel di 

coordinate n,m della prima immagine filtrata e 1' intensita 

Ing. Marco Celestino 
ABM Agenzia Br even i & Marchi 
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di tutti i pixel contenuti in un intorno del pixel di 
coordinate n,m di detta immagine di partenza, o di una 
seconda immagine filtrata da essa derivata. 

Pref eribilmente, la sommatoria delle differenze e 
5 calcolata tra la prima immagine filtrata ed una seconda 
immagine filtrata, derivata dall' immagine di partenza, 
ottenuta per ogni n,m dalla media locale calcolata su un 
intorno del pixel di coordinate n,m dell' immagine di 
partenza. 

10 Vantaggiosamente, la sommatoria delle differenze 

viene sdoppiata calcolando una sommatoria delle differenze 
positive, o deviazione positiva, e una sommatoria delle 
differenze negative, o deviazione negativa. 

Secondo 1' invenzione, il calcolo del mome 

15 centrale assoluto generalizzato dei livelli di grigio 
pixel di una immagine comprende le fasi di : 

- definizione di quattro domini circolari ©i, ©2, ©3 e 
©4, rispettivamente di raggio ri, r2, rs ed r4; 

- associazione a ogni dominio ©i di una funzione peso wi; 
20 - calcolo del valor medio pesato ^ii della mappa a livelli 

di grigio per i domini ©1 e ©2; 

- calcolo del momento assoluto centrale generalizzato 
e(n,m) con la funzione peso W3 sul dominio ©3; ; 

- scomposizione del momento centrale assoluto 

25 generalizzato e(n,m) in due filtri complementari : una!v '^^^'^%.\ 

Ing. Marco Celestino 1 '^t ]<ifA:<'^'} * 

ABM Agenzia Brex^ni <Sc Marchi "^'^^f . ' . / 
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deviazione positiva ep(n,m) e una deviazione negativa 
en(n,m) che producono, in prossimita di una discontinuita, 
due profili a campana parzialmente sovrapposti; 
- convoluzione delle due deviazioni positiva ep(n,m) e 
5 negativa en(n,m) con la funzione peso W4 sul dominio ©4; 

Uno dei principali vantaggi dell' invenzione e che 
en(n,m) produce una campana con il picco sul bordo scuro 
della discontinuita mentre ep{n,ni) produce una campana con 
il picco sul bordo luminoso. In generale, quindi^ la 
10 deviazione negativa evidenzia strutture scure su fondo 
chiaro, mentre la deviazione positiva evidenzia strutture 



chiare su fondo scuro. Inoltre, una delle due deviazioni 
evidenzia il bordo interne di una struttura mentre I'altra 
deviazione ne evidenzia il bordo esterno. Questo pud 
15 essere utile ad esempio per distinguere, una volta 
individuati, il bordo esterno da quelle interne di un vaso 
sanguigno, di una parete cardiaca, ecc. 

Un altro principale vantaggio e che la funzione 
Min (ep (n, m) ^ I en (n, m) I ) consente di superare le difficolta 
20 computazionali di ricerca degli ^^zero-crossing" della 
risposta del filtro. La zona di sovrapposizione delle due 
componenti rappresenta infatti una cresta molto stretta 
che localizza direttamente lo ^'zero-crossing", e quindi la 
discontinuita, producendo al tempo stesso un valore 
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proporzionale alia pendenza dello ^"zero-crossing". Questo 
alleggerimento computazionale rende 11 metodo vantaggioso 
per I'estrazione di contorni di immagini video organizzate 
come sequenza di fotogrammi. 



standard molto usati in computer vision come il DoG 
(differenza di gaussiane) e il GoG (gradiente di 
gaussiana) possono essere sostituiti il primo con la somma 
delle due deviazioni positiva ep(n,m) + en(n,m) e il 

10 secondo con la differenza ep(n,m) - en(n,m). Questo perche 
i risultati sono analoghi. 

Un altro principale vantaggio e che il momento 
centraie assoluto e(n,m) produce un profile a campana in 
corrispondenza di discontinuita e produce un massimo 

15 locale in corrispondenza di spigoli, di giunzioni e dei 
punti terminal! delle linee, Quindi, le prestazioni del 
momento centraie assoluto sono migliori di quelle del 
gradiente di Gaussiana proprio nei punti chiave 
dell' immagine . 

20 Un altro principale vantaggio e che il momento 

centraie assoluto e(n,m) produce un profile a campana 
anche in corrispondenza di linee e puo quindi essere usato 
per il recupero di tracciati (per esempio tracciati 
elettrocardiograf ici) da supporto cartaceo e per la loro 



5 



Un altro principale vantaggio e che due filtri 



25 conversione in segnali digitali. 
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essere acquisito con uno scanner e I'immagine ottenuta puo 
essere filtrata con il momento centrale assoluto e(n,m) in 
modo da evidenziare la traccia elettrocardiograf ica 
tramite un profile a campana. La sommita di tale profile a 
5 campana e il segnale digitale cercato e puo essere 
determinata con un algoritmo standard di localizzazione 
dei massimi locali; 

Un altro principale vantaggio e che due operazioni di 
filtraggio ottenute con il momento centrale assoluto 
10 e{n,m) ma con differenti parametri ri, r2, ra ed possono 
essere usati per implementare un algoritmo di sogliatura 
locale che mette a confront© due immagini, una con 

Ilia. V-j 's^ J- 'W -I- ^ J. .-iili'w^-i. — lll^li'-' J_ J_li.„V^il^J.^*^ **V-,-. ^_l.WV.^*4-*. v., 

minore rumore e maggiore intensita nei picchi. 
15 Vantaggiosamente, le funzioni peso sono scelte tra: 

— funzioni costanti, in cui il valore della costante e 
scelto in modo da normalizzare la funzione sul dominio ©i; 

- funzioni Gaussiane, in cui dette funzioni peso wi sono 
approssimate con funzioni normalizzate su domini circolari 

20 ©i di raggio ri=3ai. 

Secondo un altro aspetto dell' invenzione, una 
apparecchiatura per I'estrazione del contorno di strutture 
presenti in una immagine comprende almeno due filtri e una 
unita aritmetico logica (ALU) , in cui detti filtri 

25 implementano il calcolo della deviazione positiva ep (n, m^N^^j:.:. \^ 

fng. Marco Ceiestino / / (^ V^ 
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e deila deviaziune necjauiva cnvi^w^") uel riiorricntio cent-rale 

assoluto e(n,m), come sopra definite. 

In particolare, ogni filtro e caratterizzato da un 

prefissato set di parametri costituiti dai coefficienti 
5 delle funzioni peso. Tali coefficienti possono essere 

impostati in fase di inizializzazione del sistema. 

Tipicamente, il numero di bit su cui sono rappresentate le 

funzioni peso e pari al numero di bit su cui e 

rappresentata I'immagine. La ALU realizza operazioni 
10 aritmetiche e logiche sui risultati dei precedenti filtri 

in base all' operazione di filtraggio prescelta. 
Breve descrizione dei disegni 

Ulteriori caratteristiche e vantaggi dei metod 
secondo la presente invenzione, risulteranno pi 
15 chiaramente dalla descrizione che segue, fatta a titol 
esemplif icativo e non limitative, con riferimento a 
disegni annessi in cui: 

- la figura 1 costituisce una rappresentazione grafica 
dei domini 0i, 02 e 03 utilizzati nel calcolo del momento 

20 centrale assoluto del primo ordine ad un pixel 
dell' immagine; 

- la figura 2 mostra la risposta (2b) di due componenti 
Bpin.m) ed en(n,m) del primo momento centrale assoluto in 
prossimita di una discontinuita a gradino (2a) ; 

25 - la figura 3 illustra come le risposte (3a) di ep(n^m) 

Ing. Marco Celestino 
ABM A^enzia Brevetti & Marchi 
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ed en(n^m) possono essere utilmente combinate per ottenere 
different! risposte (3b-3d) ; 

- la figura 4 mostra, nel caso di due immagini di partenza 
{4a e 4d) , la risposta del momento centrale (4c e 4f) nella 

5 combinazione di figura 3b, confrontata con la risposta del 
gradiente di Gaussiana (4b e 4e) ; 

- la figura 5 mostra, nel caso di quattro iiranagini di 
partenza (5a) la risposta del filtro Min (deviazione 
positiva, I deviazione negativa | ) (5c) nella combinazione di 

10 figura 3d , confrontata con la mappa degli zero-crossing 
ottenuti come risposta ad una Differenza di due Gaussiane 
(DoG) (5b); 

- la figura 6 illustra una procedura di sogliatura locale 
secondo I'invenzione nei suoi vari stadi 6a-6g; 

15 - la figura 7 illustra una tipica applicazione della 
procedura di sogliatura locale (7d) , partendo dalle stesse 
immagini di figura 5 (7a) ; 

- la figura 8 illustra lo schema a blocchi di un apparato 
per 1' implementazione della classe di filtri proposta in 

20 grado di elaborare segnali video in tempo reale; 

- la figura 9 illustra lo schema a blocchi dei filtri di 
figura 8; 

- nella figura 10 e schematizzata 1' architettura della 
unita elementare dei quattro convolutori di fig. 9; 

25 - la figura 11 illustra 1' implementazione di un filtro 

Ing, Marco Celestino 
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bidimensionale tramite la cascata di due filtri 

monodimensionali ^^ID-FIR" ; 

- nella figura 12 e schematizzata 1' architettura 
dell'unita elementare dell' integratore di figura 9. 

Descrizione del metodo 

II metodo secondo I'invenzione parte da una 
generalizzazione del momento centrale assoluto del I ordine. 



- f(n,m) la mappa dei livelli di grigio di una immagine, 
dove f assume uno tra i 256 valori tra 0 e 255 di livelli 
di grigio per ogni pixel di coordinate n, m, 

- wi(n^m)^ W2(n,m) , W3(n,m) e W4(n,m) quattro funzioni 



ciascun punto dei domini avendo coordinate (k,l), e i 
domini essehdo scelti in funzione dell' operazione da 
realizzare . 

- rx, r2t e i rispettivi raggi di ©i, ©2/ ©3 e ©4, 

il momento centrale assoluto generalizzato del primo 
ordine e calcolato come segue: 



Indicahdo con: 




CO domini circolari ©1, ©2,- ©3 e ©4, 




(1) 



dove 



- ® e I'operatore di convoluzione. 



il valore medio /Ji(n^in) e calcolato come 




(*./)€e, 

i domini ©i sono definiti come 
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con I che rappresenta I'insieme dei numeri interi. 

Con riferimento alia figura 1, 11 primo momento 
centrale assoluto generalizzato misura la variabilita del 
valore medio dei livelli di grigio contenuti nel dominio 
circolare ©3 rispetto al valore medio locale dei livelli 
di grigio calcolato nel dominio centrale ©i; II valore 
medio dei livelli di grigio di ogni punto del dominio 
circolare ©3 viene calcolato su un intorno circolare ©2. 

Per quanto riguarda le funzioni peso, possono essere: 
— una funzione peso costante, scelta in modo da 
normalizzare la funzione sul dominio ©i. 



0 (A:,7)€0, 



con X Y^w^kJ)^! 
(*./)se/ 



(4) 



(5) 



Una funzione peso Gaussiana, approssimata con 
funzioni normalizzate su domini circolari ©i di raggio 



ri=3ai, calcolata come segue 



(kj)eei 



0 



(^,/)e0, 



(6.) 



(A:,/)«0, 



II momento centrale assoluto generalizzato e(n,m) 
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viene quindi scomposto in due filtri complementari : una 
deviazione positiva ep(n,m) e una deviazione negativa 
enin^m) , e precisamente : 



{n. m) = {n, /w) ® J] m)- ^^{n-k.m- l))w^ {k, /) 



3P 



(7) 



e„ («, m) = («, m) ® 5] J] (//^ (w, m)- /d^in- k,m- l))w^ (k, /) 

5 dove i domini ©3p e 03n sono definiti come: 
®3p ={(^,/)e©3 :/i,(n,;n)>/i2(«-^,/n-/)} 

(8) 

©3n ={(A:,/)€©3 :^i,{n,m)<fi^{n-k,m-l)} 

II momento centrale assoluto generalizzato puo 
quindi essere ottenuto come: 

10 e(n,m) = ep(n,m) - e„(n,m) (9) 

La figura 2 mostra la risposta (2b) di en(n^m) 
ep(n^m) ad una immagine test (2a) con una discontinuita\^ 
ideale a gradino tra due zone a diverso livello di grigio 
20 e 30. In corrispondenza della discontinuita 10, en('n,/n; 

15 ed ep(n^m) producono due profili a campana parzialmente 
sovrapposti., Mentre en(n^m) produce una campana con il 
picco sul bordo scuro della discontinuita, ep(n^in) produce 
una campana con il picco sul bordo luminoso. In generale 
la deviazione negativa evidenzia strutture scure su fondo 

20 chiaro, mentre la deviazione positiva evidenzia strutture 
chiare su fondo scuro. 

'La figura 3 mostra schematicamente come le uscite 




ing. Marco Celestino 
ABM Anemia Brevetti <&. Marchi 
Isctii»ck<kl'dib6\/. 544 



PI POO? A 0 0 0 0 5 1 



- 18 - 



dei filtri epfn^m) ed en(n^m) possono essere combinate per 
ottenere alia discontinuita different! risposte. 
- Fig- 3a: le singole componenti ep(n^m) e \en(nrm)\ sono 
mostrate come curve a campana separate. 
5 - Fig. 3b: La somma delle componenti positiva e negativa 
ep(n,m) + \en(n^m)\ da un singolo profilo a campana il cui 
picco localizza la discontinuita. Questa risposta e molto 
simile a quella che potremmo ottenere con un gradiente di 
Gaussiana . 

10 - Fig- 3c: La differenza delle componenti positive e 
negative ep(n^m) - \en(n,m)\, invece, produce una risposta 
simile a quella che potremmo ottenere con un filtro 
Lapiaciano di Gaussiana- II filtro produce infatti uno 
zero-crossing in corrispondenza della discontinuita. 

15 - Fig. 3d: Le difficolta computazionali di ricerca degli 
zero-crossing della risposta del filtro vengono superate 
con la funzione Min(ep(n,m) ^\en(n^m) \ ) , La zona di 
sovrapposizione delle due componenti rappresenta infatti 
una cresta molto stretta che localizza lo zero-crossing, e 

20 quindi la discontinuita, producendo al tempo stesso un 
valore proporzionale alia pendenza dello zero-crossing. 

Generazione di un massimo locale alle discontinuita 

Con la somma ep(n,m) -h \en(nrm)\ si ottiene la 
funzione e(n^m)^ la quale produce un valore nullo sopra le 



25 regioni omogenee ed un massimo locale in corrisponden 
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delle discontinuita, II momento centrale assoluto e un 
indice di dispersione che produce una mappa simile a 
quella prodotta dal gradiente di Gaussiana. La figura 4 
mostra come, a partire da due immagini (4a e 4d) la 

5 combinazione lineare ep(n,m) -h \en(n,m)\ genera (4c e 4f) 
un picco in corrispondenza della discontinuita . tuttavia 
importante osservare che, diversamente dal gradiente di 
Gaussiana (4b e 4e) , il momento centrale assoluto puo 
produrre un profile a campana sia in corrispondenza di 

10 linee che in corrispondenza di discontinuita a gradino e 
che puo produrre un massimo locale in corrispondenza di 
spigoli, di giunzioni e dei punti terminali delle linee • 
Quindi, le prestazioni del momento centrale assoluto sono 
migliori di quelle dei gradiente di Gaussiana proprio nei 

15 punti chiave dell' immagine . 

Nel caso di funzioni peso Gaussiane, le risposte del 
momento centrale assoluto risultano diverse al variare 
della dimensione dei domini ©i, ©2, ©3- In particolare il 
momento centrale assoluto generalizzato produce, 

20 rispettivamente, in caso di : 

- cri<a2<(y3 - un solo profile a campana alle discontinuita; 

- (Ji<a2<l/2(j3 massimi locali agli spigoli; 

- <Ti<b2<2/3crj ; 6o'2>d ; o'3>d ~ una campana a linee di 
ampiezza d. 

25 L' introduzione della funzione peso wi(n,m) impedisper^. - ' 

•'».•.' 
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perd all' operatore di evidenziare tutte le giunzioni con 
un massimo locale. D' altra parte la funzione peso W2(n^m) 
e necessaria per contrastare gli effetti del rumore 
normalmente presente nelle immagini e non puo essere 
eliminata. Per evidenziare tutte le giunzioni con un 
massimo locale dobbiamo mantenere separate le due 
deviazioni positiva ep(n^m) e negativa en(nrm) in quanto 
esse forniscono un massimo locale a tutte le giunzioni. 

Vale la pena notare che il momento centrale di 
ordine uno nella sua forma standard e identicamente uguale 
a zero. Nel nostro caso, tuttavia, I'operatore ottenuto e 
una generalizzazione del momento centrale di ordine uno e 
non e identicamente uguale a zero. II momento centrale di 
ordine uno: 



e equivalente all'eq.(l) senza valore assoluto e 

come tale puo quindi essere ottenuto come somma delle due 

componenti negativa e positiva dell'eqd). Sviluppando 

I'eq. (10) si ottiene: 
c{n, m) = (/7, w) + e„ («, m) = 

= / m) (8) (w^ («, m) 0 (n, m) - («, m) <S> (n, m) ® {n, m)) = (11) 

= m) ® (w^ (n, m) - {n, m)) 

L'eq.(ll) mostra come nel caso di funzioni peso 
Gaussiane la somma ep(nrm) + en(nrm) fornisca una risposta 
equivalente a quella fornita da un filtro DoG (differenza 



(10) 
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di Gaussiane) • 

Mappa e intensita degli zero-crossing 

Qui di seguito viene indicato un filtro ottenuto da 
un semplice algoritmo per ottenere contemporaneamente la 
mappa degli zero-crossing e una stima della loro intensita 
II momento centrale assoluto pud essere diviso nelle 
componenti positiva e negativa, e quattro funzioni peso 
possono essere introdotte nei filtri originali. Nel caso 
di funzioni peso Gaussiane la somma delle risposte delle 
due componenti positiva e negativa e equivalents alia 
risposta di un filtro DoG, indipendentemente dall 
aperture Oi, 02, Oj, a^. Tuttavia, nel punto in cui le d 
componenti della risposta del DoG attraversano lo zero, le 
due creste generate dalla deviazione negativa e positiva 
del momento assoluto si sovrappongono parzialmente . II 
profilo dell' area di sovrapposizione e una sottile cresta 
il cui picco localizza la discontinuita . L'altezza del 
picco inoltre e proporzionale al valore di contrasto 
rappresentato dalla discontinuita. Quindi una semplice 
funzione algebrica, quale la funzione Min (deviazione 
positiva^ I deviazione negatival), puo fornire 

contemporaneamente una mappa Mpn simile alia mappa degli 
zero-crossing prodotta con il metodo di Marr-Hildreth 
(fig. 3) e una stima dell ' intensita degli zero-crossing. 



In fig. 5 " sono rappresentate quattro immagini (5a), 
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di cui la prima e una immagine sintetica con rumore 
gaussiano e le altre sono tre generiche immagini di 
esterni. Sono quindi confrontate le mappe (5b) degli zero- 
crossing ricavate dall' applicazione di un filtro DoG 
5 (metodo di Marr-Hildreth) e le mappe (5c) ottentute con la 
procedura Min (deviazione positiva, Ideviazione negatival) 
secondo 1' invenzione; queste ultime mostrano contorni 
simili a quelli delle mappe di zero-crossing. Nell'esempio 
le impostazioni dei filtri passa basso sono le stesse per 

10 entrambi gli approcci . Sono state usate le classiche 
aperture a2=1^6a'i 0^=0 di un filtro DoG. Una Gaussiana 
W4(n^ m) con apertura a4=0.7 pixel e stata infine usata 
per assicurare la sovraooosizione delle due deviazioni del 
momento assoluto. Se questa conf igurazione e utilizzata le 

15 due deviazioni sono uguali alle due componenti positiva e 
negativa della risposta di un filtro DoG quando sono 
filtrate separatamente con una Gassiana di apertura uguale 
a 0.7 pixel . 

Procedura locale di sogliatura 
20 Una procedura locale di sogliatura puo essere 

ottenuta combinando due mappe prodotte dall'uso dei 
seguenti filtri: 



e"(Ai,/w) = w,(n,m)®2 Y}firi,m)-f{n'k,m--l)\w,{kJ) (13) 



(*,/)€e3 
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In presenza di rumore, quando rj = r^, su zone 
omogenee dell ' immagine I'equazione (12) produce un livello 
di rumore inferiore a quello prodotto dall' eq. (13) . In 
corrispondenza dalle discontinuita si ottiene invece il 
5 risultato opposto. Alle discontinuita I'eq. (12) produce 
profili a campana piu alti di quelli prodotti 
dall'eq. (13) . Dato che I'eq. (12) produce massimi locali 
piu alti alle discontinuita e valori di rumore piu bassi 
sulle zone omogenee rispetto a quanto prodotto 

10 dall' eq. (13) ^ quest'ultima puo essere utilizzata come 
soglia della risposta ottenuta con 1' eq. ( 12) . 

Con riferimento alia figura 6, quindi, partendo da 
un' immagine test 6a, alia, quale e stato sommato rumore 
gaussiano, vengono applicati i filtri relativi alle eq. 

15 (12) e (13), ottenendo rispettivamente le immagini 6b e 
6c; i grafici 6d ed 6e mostrano.il profile dei livelli di 
grigio ottenuti rispettivamente dalle immagini 6b e 6c in 
corrispondenza delle righe qui rappresentate; sogliando il 
grafico 6d con la curva del grafico 6e si ottiene il 

20 profile 6g; in • quest ' ultimo e evidente 1' eliminazione del 
rumore quale effetto della procedura di sogliatura locale. 

Come mostrato in figura 7, applicando la procedura 
di sogliatura , locale alle immagini di riga 7a, le stesse 
di figura 5a, filtrandole con I'eq. (10) si ottengono le 



25 .immagini della riga 7b, filtrandole con I'eq. (11) si<^ 
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ottengono le iinmagini della riga 7c, e sogliando le 
immagini della riga 7b con le immagini della riga 7c si 
ottengono le immagini della riga 7d. 

Da quanto sopra, e quindi evidente che i filtri 
secondo 1' invenzione, derivati dalla generalizzazione del 
momento centrale assoluto del primo ordine, pur presentando 
analogie funzionali con alcuni filtri correntemente 
utilizzati (gradiente di Gaussiana, Laplaciano di Gaussiana, 
differenza di Gaussiane, ecc), sono vantaggiosi perche in 
un unico filtro si trovano concentrate le seguenti 
proprieta: 

- al variare del parametri del momento centrale assoluto 
del prirrio ordine i filtri oLLenuti possono fornire gli 
stessi risultati dei filtri tradizionali quali gradiente 
di Gaussiana, Laplaciano di Gaussiana e differenza di 
Gaussiane, 

- si possono implementare filtri derivativi sia del primo 
e del secondo ordine, 

- possono essere usati per ottenere al tempo stesso la 
localizzazione degli zero-crossing della derivata seconda 
e una stima della loro importanza, 

possono essere usati per 1' implementazione di un 
algoritmo di sogliatura locale, 

-■ e possibile evidenziare separatamente il bordo interne 

e quello esterno di una struttura, o separatamente ogqettx^-^^rZ^r- 
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scuri su fondo chiaro e oggetti chiari su fondo scuro, 

- e possibile evidenziare con massimi locali isolati 
punti chiave dell' immagine quail spigoli, giunzioni e 
punti terminali di linee o punti isolati, 

5 - e possibile evidenziare la linea centrale di linee 
aventi spessore non nullo con una cresta, e essere usati 
per estrarre tracciati (per esempio ECG) da supporti 
cartacei, 

- nonche possono essere usati, tramite i suddetti punti 

10 chiave, per ricostruire immagini panoramiche di una 

struttura da una pluralita di scansioni prese da diver^^^ ' 

I ^ <^ 

punti di vista. Y^*^ 
o Q V- j_ ^ T_ QT-^ g Ht_ ijrj H zi p htaizx^^pp^ oreferito ^ 

Nelle applicazioni dove e richiesto il filtraggio di 

15 poche immagini, I'algoritmo proposto puo essere 

implementato ef f icientemente su un comune calcolatore. In 

applicazioni piu onerose, quali il filtraggio in tempo 

reale di segnali video, la classe di filtri puo essere 

implementata tramite un dispositive hardware dedicate come 

20 quelle decritto in quest ' esempio . 

L' apparato proposto, in base alia presente 

invenzione, e un dispositive conf igurabile in grade di 

implementare 1' intera classe di filtri. II sistema e in 

grado di elaborare in tempo reale un flusso continue 

25 d' immagini con risoluziene e frame-rate tipici del comuni 
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standard video. Sia il segnale video di ingresso sia 
quello di uscita sono in formato digitale.- Ad esempio, 
essi possono essere una sequenza d'interi ad 8 bit 
rappresentanti immagini su 256 livelli di grigio. Nel caso 
in cui la dimensione delle immagini sia di 512x512 pixels 
e si abbiano 25 fotogrammi al secondo, allora la frequenza 
di pixel, cioe la velocita dei dati, e 
fp^5 12x512x25=6. 5536 MHz. 

L' apparato pud essere descritto tramite un 
linguaggio HDL (Hardware Description Language) che ne 
definisce in modo univoco 1' architettura e ne permette 
1' implementazione su di un circuito integrate. La scelta 
del dispositivo su cui mappare 1' architerrura dipende, 
invece, dalle esigenze poste in fase di realizzazione . Nel 
caso si vogliano privilegiate flessibilita e rapidita di 
sviluppo del prodotto e possibile utilizzare dispositivi 
di logica programmabile, quali FPGA (Field Programmable 
Gate Array) e CPLD (Complex Programmable Logic Device) . Se 
viceversa si intendono privilegiare le prestazioni ed il 
costo finale del dispositivo, la soluzione da adottare pud 
essere quella di un circuito integrate ASIC (Application- 
Specific Integrated Circuit) . In entrambi i casi puo 
essere indispensabile I'utilizzo di banchi di memoria 
esterni, dove memorizzare temporaneamente i dati. 

L' architettura dell' apparato, illustrata in figu 

fng: Marco Celestino 
ABM Agenzia Brevetti & Marchi 



PI 9009 A 0 0 0 0 5 1 



8, e costituita da uno o piu filtri 1 e da una ALU (unita 
aritmetico logica) 2. I filtri implementano il calcolo 
della deviazione positiva ep(n,m) e della deviazione 
negativa en^n^in;. Ogni filtro e caratterizzato da un 

5 prefissato set di parametri costituiti dai coefficienti 
delle funzioni peso wx(k^l) W2(k,l) wsik^l) W4(krl). Tali 
coefficienti possono essere impostati in fase di 
inizializzazione del sistema. Tipicamente, il numero di 
bit su cui sono rappresentate le funzioni peso e pari al 

10 numero di bit su cui e rappresentata I'immagine. La ALU 
realizza operazioni aritmetiche e logiche sui risultati 
dei precedenti filtri in base all' operazione di filtraggio 
prescelta . 

La conf igurazione minima dell' apparato e quella in 
15 cui si ha un unico filtro e dove la ALU realizza una 
semplice operazione fra deviazione positiva e deviazione 
negativa. Nel caso, ad esempio, in cui si voglia ottenere 
un massimo locale in corrispondenza delle discontinuitav 
1' operazione effettuata dalla ALU e la differenza fra le 
20 due deviazioni. 

g(«, m) = m) = («, m) -e^in.m) (14) 
Dalla somma delle due deviazioni: 

g{n, m) = e^in,m) + e„{n,m) (15) 
si ottengono, invece, zero-crossing in corrispondenza 
25 delle discontinuita. 
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Nel caso piu generale in cui. si abbiano due filtri, 
la ALU ha in ingresso le componenti ottenute tramite il 
primo filtro e^p(n^m) e^ni^fTa) e quelle ottenute tramite 
il secondo filtro e''^p(n^m) e^^n(^/^) • Tale conf igurazione 
pud essere usata, ad esempio, per operazioni di sogliatura 
locale, L'operazione effettuata dalla ALU e la seguente: 

\ (16) 



, . . V. , e\n,m)>e'\n,m) 
g(n,m) = \^ (17) 



L' architettura del filtro 1, illustrata in figura 9, 
e costituita da quattro convolutori FIR bidimensionali e 
da un integratore. I quattro convolutori imp 1 erne ncano le 
convoluzioni (2) dell' immagine f(n^m) con le maschere 
wi(krl) W2(k,l) e le due convoluzioni (7) dei segnali 
ep(nrm) ed en(n^m) con la maschera W4(k,l). 

L' integratore implementa le sommatorie contenute 
nell' equazione (7) e riportate di seguito: 



_ , (18) 

Nel caso in cui si elaborino iinitiagini su 256 livelli 
di grigio, e possibile rappresentare tutti i segnali 
f(nrm)j. ]ji(n,m) , iJ2(nrm), ep(n^m) , en(n,m) , ep(n,m) ^ 
en(n^m)^ tramite numeri interi ad 8 bit. Questa 
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approssimazione costituisce un buon compromesso* fra 
precisione del risultato e complessita dell' hardware . 

I convolutori sono costituiti da un certo numero 
d'unita elementari 11 cui schema a blocchi e illustrato in 
5 figura 10. L'operazione implementata da ciascuna unita e 
la moltiplicazione-e-somma, owero: 

- moltiplicazione fra I'elemento f(n^m) dell' inunagine ed 
11 coefficiente w(k^l) della maschera di convoluzione . 

- somma s(krl) fra 11 risultato della moltiplicazione ed 
10 11 valore s(k,l-l) ottenuto dall' unita precedente. > 



- un registro 6 dove e memorizzato 11 coefficiente della 
maschera di convoluzione. 

15 - un moltiplicatore 7. 

- un sommatore 8 . 

Nel caso in cui si elaborino immagini su 256 livelli 
di grigio, gli element! 5 e 6 sono implementati tramite 
registri ad 8 bit. Sia il moltiplicatore 7 che 11 

20 sommatore 8 sono invece implementati utilizzando logica 
combinatoria. II primo ha ingressi ad 8 bit ed uscita a 16" 
bit mentre la dimensione del secondo dipende dalla 
dimensione della somma parziale s(k^l-l). In alcune 
implementazioni pud essere vantaggioso realizzare 

25 moltiplicatore e sommatore tramite un'unica rete 



un'unita di ritardo 5. -BPg^ 



Gli element! di cui e costituita ciascuna unita sono/'^^^ 
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combinatoria avente tre ingressi ed un'uscita. 

Se la velocita d' elaborazione dell'unita e pari alia 
frequenza di pixel, e possibile realizzare una 
moltiplicazione-e-sonutia per ogni ciclo di clock. Questo 
significa che per realizzare la convoluzione in tempo 
reale e necessario utilizzare un numero di unita pari alia 
dimensione della maschera di convoluzione w(kfl) . 
Tuttavia, se si utilizzano dispositivi a semiconduttore 
abbastanza veloci, e possibile ridurre notevolmente la 
dimensione del convolutore e quindi la complessita del 
circuito integrato che lo implementa. Si consideri il caso 
in cui la velocita d' elaborazione dell'unita sia N volte 
j-ci ^ i ^•4w;;:=ii<ici j. i i i^ci-L caso G possiDile utilizzare - 

la stessa unita per realizzare la moltiplicazione-e-somma 
fra I'immagine ed N valori della maschera di convoluzione 
w(k,l), i quali saranno caricati ciclicamente nel registro 
6. Utilizzando quest ' accorgimento, la dimensione del 
convolutore si riduce di un fattore N, 

Tramite un convolutore bidimensionale e possibile 
implementare la convoluzione con una qualsiasi funzione 
peso Wi(k,l) definita nell'intorno circolare ©i. Tuttavia 
nel caso si utilizzino funzioni particolari, come ad 
esempio le Gaussiane, e possibile introdurre un secondo 
accorgimento per ridurre la dimensione del convolutore. 

Si considerino funzioni peso Wi(k,l) definite su 
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dominio quadrate: 

0/ = {(^,/)e/2:-r,<^<r,;-r, </</;} (19) 

Dal confronto con la (3) si puo facilmente dedurre 
che il dominio quadrato include quelle circolare : 0^ e 0/ . 
Nell' ipotesi in cui i valori numerici della funzione 
Wi(k^l), siano trascurabili al di fuori dell' intersezione 
fra i due domini, la convoluzione con Wi(k,l) e la 
convoluzione con Wi(kfl) producono risultati equivalent i. 

Se, inoltre, la funzione Wi(k^l) e separabile, cioe 
e possibile definire due funzioni Wii(k) e w±2(l) tali che: 

M^Xk^)^w,,{kyw,,{l) (20) 

allora la convoluzione bidimensionale con Wi(k^l) e 

equivalente alia serie di due convoluzioni 

monodimensionali con Wii(k) e Wi2(l) rispettivamente . Nel 

caso, ad esempio, del calcolo dei valori medi della 

funzione f(n^m), I'equazione (2) diventa: 

( \ 

con Wii(k) e Wi2(l) definite sui domini monodimensionali: 
M II Hi ^22) 

Con queste ipotesi, il convolutore bidimensionale 
pud essere implementato dalla cascata di due convolutori 
monodimensionali, come illustrate in figura 11. Nel caso 
in cui i dati siano ad 8 bit, la dimensione degli elementi 
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dell'unita elementare del primo convolutore 9.1 rimane 
quella descritta precedentemente . I dati r(n^m) in uscita 
dal primo convolutore sono, ihvece^ rappresentati su un 
numero maggiore di bit: 16 Log2 (2 - rn-hl) . Tuttavia, nella 
maggior parte dei casi presi in esame, e possibile 
approssimare tale dato tramite un intero ad 8 bit. Con 
quest ' ulteriore approssimazione anche la dimensione degli 
elementi dell'unita elementare del secondo convolutore 9.2 
rimane invariata. 

Nel passaggio da una convoluzione bidimensionale con 
dominie quadrato a due convoluzioni monodimensionali , il 
numero di unita elementari necessarie all' operazione si 
riQuce da (^ri~hl) - a ^(2ri-hl). Nel caso di domini 
circolari, il numero di unita elementari e generalmente 
minore di (2ri-hl)^ essendo c©i. Tuttavia, il vantaggio 
che si ottiene nel passaggio a convoluzioni 
monodimensionali e ancora notevole. 

L' integratore e costituito da unita elementari la 
cui architettura e riportata in figura 12. Con riferimento 
alle equazioni (18), le operazioni implementate 'da 
ciascuna unita sono: 

- sottrazione fra i valori }2i(n^m) ed 112 (n-k .m-l) . 

- moltiplicazione con il coefficiente della matrice 

W3(k,l) . 

- se il risultato della sottrazione e positive: somma 
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Sp(k,l) del risultato della moltiplicazione con il valore 
Sp(kfl-l) ottenuto dall'unita precedente. 

- se il risultato della sottrazione e negativo: somma 
Sn(k,l) del risultato della moltiplicazione con il valore 

5 Sn(k,l-1) ottenuto dall'unita precedente. 

Gli elementi di cui e costituita ciascuna unita 
elementare sono : 

- un' unita di ritardo 10. 

- un registro 11 dove e memorizzato il coefficiente della 
10 maschera di convoluzione . 

- un sottrattore 12 con riporto (borrow) . 

- un moltiplicatore 13. 

- una costante zero 14. 

- due multiplexer 15.1 e 15.2. 
15 - due sommatori 16.1 e 16.2. 

Nel caso in cui si elaborino immagini su 256 livelli 
di grigio, gli elementi 10 e 11 sono implementati tramite 
registri ad 8 bit. Tutti gli altri elementi, numerati da 
12 a 16.x, sono invece implementati tramite logica 
20 combinatoria . La sottrazione 12 e fra operandi ad 8 bit 
con risultato su 8 bit mentre la moltiplicazione 13 e fra 
operandi ad 8 bit con risultato su 16 bit. La dimensione 
dei multiplexer 15.x e dei sommatori 16.x dipende dalla 
dimensione delle somme parziali Sp(krl-l) ed Sn(k^l-l) . 

25 - Se la . velocita d' elaborazione dell' unita e pari alia 
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frequenza di pixel, 1' integratore e costituito da Af^ unita 
elementari, dove M e la dimensione della maschera W3(k^l). 
Analogamente a quanto affermato per i convolutori, 
disponendo di un dispositive a semiconduttore veloce e 
5 possibile utilizzare la stessa unita elementare per 
calcolare N elementi della sommatoria. In tal mode anche 
la dimensione dell' integratore si riduce di un fattore N. 

La descrizione di cui sopra di una forma esecutiva 
specif ica e in grado di mostrare I'invenzione dal punto di 

10 vista concettuale in modo che altri, utilizzando la 
tecnica nota, potranno modificare e/o adattare in varie 
applicazioni tale forma esecutiva specifica senza 
ulterior! ricerche e senza allontanarsi dal concetto 
inventive, e, quindi, si intende che tali adattamenti e 

15 modifiche saranno considerabili come equivalent! della 
forma esecutiva esemplif icata . I mezzi e i material! per 
realizzare le varie funzion! descritte potranno essere di 
varia natura senza per questo uscire dall'ambito 
dell' invenzione . Si intende che le espressioni o la 

20 terminologia utilizzate hanno scopo puramente descrittivo 
e per questo non limitative 
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1. Un metodo automatico di localizzazione ed estrazione 
del contorno di strutture present! in una immagine di 
partenza, caratterizzato dal fatto di prevedere il 
filtraggio di detta immagine di partenza tramite il 
momento centrale assoluto e(n^m) dell' intensita dei pixel 
di detta immagine, essendo n e m le coordinate di ciascun 
pixel, in cui detto momento centrale assoluto e ottenuto 
con le seguenti fasi: 

- determinazione per ogni n^m della media locale 
calcolata su un intorno del pixel di coordinate n^m 
dell' immagine di partenza, ottenendo una prima immagine 
f iltrata; 

determinazione per ogni n,/n della sommatoria deiie 
differenze in valore assoluto tra 1' intensita del pixel di 
coordinate n,m della prima immagine filtrata e 1' intensita 
di tutti i pixel contenuti in un intorno del pixel di 
coordinate n,in. di detta immagine di partenza, o di una 
seconda immagine filtrata derivata da detta immagine di 
partenza . 

2. iVIetodo secondo la rivendicazione 1, in cui detta 
sommatoria delle differenze e calcolata tra detta prima 
immagine filtrata e detta seconda immagine filtrata, in 
cui detta seconda immagine filtrata e ottenuta per ogni 
n,/n dalla media locale calcolata su un intorno del pixel 



di coordinate n^m di detta immagine di partenza. 
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3. Metodo secondo la rivendicazione 1, in cui detta 
sommatoria delle differenze viene sdoppiata calcolando una 
sommatoria delle differenze positive, o deviazione 
positiva, e una sommatoria delle differenze negative, o 
deviazione negativa . 

4. Metodo secondo la rivendicazione 1, in cui detto 
momento centrale assoluto e(n,m) e calcolato in modo 
generalizzato come segue: 

ein, m) = («, /w) ® J] X rn)- /j^{n-k,m- l)\w^ {k, I) ( 1 ) 

(*./)e03 

essendo n^m le coordinate dei pixel di una mappa f(n^m) 
di detta immagine; wi(nrm)r W2(n,m), W3(n^m) e W4(n,m) 
quattro funzioni peso definite su quattro domini circolari 
©1/ 02/ 03 e 04, rispettivamente di raggio ri, r2, rj e e 



definiti come 0^ = {(A:,/) e : VA:^ -h/^ <r.} ; ® e I'operatore di 
convoluzione; /ii(«,m) = ^ ^f{^^^)'^\{n,m) il valore medio sul 

(*,/)e0, 

dominio .0i di detta mappa e costituente detta prima 
immagine filtrata; /i2(«,/w) = ^ ^fin-k,m-l)w^{kj) il valore 
medio sul dominio 02 di detta mappa e costituente detta 
seconda immagine filtrata. 

5, Metodo secondo le rivendicazioni 3 e 4, in cui a 
partire da detto momento centrale assoluto generalizzato 
sono utilizzati come ulteriori f iltri detta deviazione 
positiva ep(n,m) e detta deviazione negativa en(n^itt) come 



/if 
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ep («, m) = (m, m) ® ^ J] («, m) - fi^{n - k,m - l))w^ {k, I) 
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(*./)e0 



= ® Zi Z C"i ('^^ m)-]u^{n-k,m- l))w^ (A, /) 

dove i domini ©3p e ©3n sono definiti come: 

03« ={(*,/) e ©3 :^i,{n,m)<fi^{n-k,m~l)} 



(7) 



(8) 



6. Metodo secondo la rivendicazione 4, in cui detto 
5 calcolo del momento centrale assoluto generalizzato 
dell' intensita del pixel prevede le fasi di : 

- definizione di detti domini circolari ©1, ©2, ©3 e ©4, 
nell'intorno di ciascun punto dell' immagine di partenza 
in cui ©1, ©3 e ©4 sono centrati su n, m e ©2 e un domini' 

10 centrato rispettivamente su ogni punto interno a ©3; 

- associazione a ciascun dominio ©i, con i che va da^ 1 a 
4, di una funzione peso Wi e calcolo del valor medio 
della mappa a livelli di grigio per i domini ©1 e ©2 sulla 
base di detti Wi; 

15 - calcolo del momento assoluto centrale generalizzato 
e(n,m) sulla base della funzione peso W3 sul dominio ©3; 

- scomposizione del momento centrale assoluto 
generalizzato efn,in>) in una deviazione positiva ep(n^m) e 
una deviazione negativa en(nrTn) producendo in 

20 corrispondenza di detta discontinuita due profili a 
campana parzialmente sovrapposti; 

- convoluzione delle due deviazioni positiva ep(nfm) 
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negative en(n,m) con la funzione peso W4 sul dominio ©4; 

- misura dell'entita della discontinuita utilizzando come 
filtri dette deviazione positiva ep(n,m) e deviazione 
negativa en(n,m). 

5 7. Metodo secondo la rivendicazidne 6^ in cui detta misura 
dell'entita della discontinuita e fatta mediante 
I'estrazione di una funzione 
Mln (Bp (n,m) , | (riria) \ ) 

8. Metodo secondo la rivendicazione 6, in cui detta misura 
10 dell'entita della discontinuita e fatta con la sottrazione: 

e{n, m) = («, m) - e„ (/?, m) 

9. Metodo secondo la rivendicazione 6, in cui un filtro 
. DoG (differenza di gaussiane) si ottiene utilizzando la 

somma di dette deviazioni positiva ep(n,m) e negativa 
15 en{n,m) del momento centrale assoluto, 

10. Metodo secondo la rivendicazione 1, in cui dette 
immagini di partenza sono: 

— immagini biomediche, ottenute con impulsi ultrasonici, 
PET, SPECT, CAT, MR, etc, tra cui immagini anatomiche, o 

20 immagini di funzione, ottenute per mezzo di sequenze 
temporali di viste anatomiche di una zona particolare di un 
organo, o immagini di perfusione, ottenute sullo stesso 
organo dopo trattamento del paziente con sostanze che 
mettano in risalto la perfusione nell' organo. ^--r 




PI 2009 AO 0 0 0 5 1 

- 39 - 

— immagini di tracciati acquisite con uno scanner con lo 
scopo di convertire il tracciato su carta in segnale 
digitals. 

ll.Metodo second© la rivendicazione 10 in cui, nel caso di 
5 immagini di tracciati, I'immagine ottenuta e filtrata con 
detto memento centrale assoluto e(n,m) evidenziando una 
traccia tramite un profilo a campana, la cui sommita e il 
segnale digitale cercato, essendo prevista una ulteriore 
fase di determinazione del segnale digitale con un 
10 algoritmo standard di localizzazione del massimi locali. 

12.iyietodo di estrazione del contorni, secondo la 
rivendicazione 4, caratterizzata dal fatto che in 
corrispondenza di una discontinuita, detto momento 
centrale assoluto generalizzato, calcolato come 
15 e\n,m) = J^ j;|/i,(«,m)-/(n-/c,m-/)|>V3(^,/) (12) 

(/:,/)603 

viene confrontato con una soglia derivata da detto 
momento centrale assoluto generalizzato, calcolata come 
(n, aw) = w, («, /w) ® X Z |/(^' -f(n-k,m- l)\w, (kj) (13) 

(A./)ee3 

13 . Apparecchiatura per 1' estrazione del contorni da 
20 immagini video organizzate come sequenza di fotogrammi, 
caratterizzata dal fatto di utilizzare una unita 
aritmetico logica ed uno o piu filtri che implementano le 
operazioni di deviazione positiva ep(n,m) e di deviazione 
negativa en(n,m) del momento centrale assoluto 
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generalizzato dell' intensita del pixel di detta imniagine 
definiti come: 

ep («, m) = (n, /w) ® ^ ^ («, m)- fi^in-k.m- l))w^ (k, /) 

(7). 

e„ («, w) = (rt, w) ® ^ ^ (/i, (w, m) - (« - ^, /n - /))w3 0 
dove i domini 03p e 03n sono definiti come: 

5 (8) 

14 . Apparecchiatura secondo la rivendicazione 13 in cui 
detti filtri sono costituiti da quattro convolutori 
bidimensionali e da un integratore. 

15 . Apparecchiatura, secondo la rivendicazione 14, in cui 
10 ognuno di detti convolutori bidimensionali e sostituito 

dalla cascata di due convolutori monodimensionali . 

p.p. CNR CONSIGLIO NAZIONALE DELLE RICERCHE; ESAOTE S.p.A. 
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